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Abstract 

By using AdS/CFT, we study the critical electric field, the Schwinger pair creation rate 
and the potential phase diagram for the quark and anti quark in four confining supergravity 
backgrounds which are the Witten QCD, the Maldacena-Nunez, the Klebanov-Tseytlin and 
the Klebanov-Strassler models. We compare the rate of phase transition in these models 
and compare it also with the conformal case. We then present the phase diagrams of the 
entanglement entropy of a strip in these geometries and find the predicted butterfly shape in 
the diagrams. We find that the phase transitions have higher rate in WQCD and KT relative 
to MN and KS. Finally we show the effect of turning on an additional magnetic field on the 
rate of pair creation by using the imaginary part of the Euler-Heisenberg effective Lagrangian. 
The results is increasing the parallel magnetic held would increase the pair creation rate and 
increasing the perpendicular magnetic held would decrease the rate. 
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1 Introduction 


Schwinger effect in quantum field theory |l] is the creation of pairs of particles in the presence of 
a strong electric or magnetic held. Intuitively when the strength of the electric held reaches the 
string tension between the quark anti-quark, it can break the string that is attaching them and 
so the virtual pairs can become on-shell and a current can be created. Semenoff and Zarembo |2| 
hrst studied the Schwinger effect holographically. Then the authors in |3| and |4], using the 
Nambu-Goto action on a probe brane and by calculating the free energy and using the hrst law 
of thermodynamic, calculated the entanglement entropy of a quark and an anti-quark which are 
accelerating in an electric held in the AdSs background. A similar calculation in a different setup 
was also done in (5]. So one might want to try this in other backgrounds and specifically in the 
confining ones. 

On the other hand, there are some ideas on how to measure the entanglement entropy in the 
lab by using the huctuations of current which is flowing through a quantum point contact as the 
probe, E 0. which still are not quite successful. Another conjecture is that there might be a 
relation between the Schwinger effect and the entanglement entropy. Finding such a relation could 
be a breakthrough since the Schwinger pair creation rate could also act as an entanglement meter 
in the condensed matter systems. To explore for such a relation, in the hrst step, we study the 
phase diagrams of Schwinger effect and of entanglement entropy for four confining geometries; the 
Witten-QCD, the Maldacena-Nunez, the Klebanov-Strassler and the Klebanov-Tseytlin which are 
dual to Af = 1 held theories. 

For calculating the entanglement entropy for accelerating particles one can use the method 
in [4], For doing so, hrst one needs to find the string world sheet prohlc in these confining back¬ 
grounds. For the case of AdS geometry, due to the large amount of symmetries, the corresponding 
PDE equation of motion is simple and has been solved in [§]. Mikhailov also found a simple 
linear relation [9] to find the string world sheet prohle based on the position of the quark and 
anti-quark on the boundary that can only be used for the AdSs background. For other geometries 
such as the confining backgrounds there is no Mikhailov-like equation and for finding the string 
prohle one needs to solve several more difficult PDE equations analytically which for our super¬ 
gravity background geometries we present them in section 2. We just do the similar calculation 
for the Minkowski background and we find the free energy of the quarks and anti quarks in a hat 
background in section 3. 

I 11 order to look for a relationship between the phase transitions of Schwinger effect and the 
entanglement entropy, we study the phase transitions in an electric potential similar to the pro¬ 


cedures in 11 10 12 13 14 15 in section 4. By using DBI action, we calculate the critical 


electric helds in these geometries where the potential becomes catastrophically unstable and we 
then find the phase diagrams numerically. Interestingly the phase diagrams of all of these geome¬ 
tries are very similar where three different phases can be detected as has been predicted in 10 as 
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a universal feature of all confining geometries. We also compare the diagrams with the conformal 
case of Klebanov-Witten which only two phases can be detected. 

In section 5 we look at the phase diagrams of the entanglement entropy of a strip, similar to the 
and (181. Klebanov, Kutasov and Murugan found a generalization of Ryu- 


calculation in 16 


Takayanagi relation for the non-conformal geometries 16 . Then in [17], the authors presented the 
plots of the phase diagram for geometries constructed by Dp brane compactified on a circle which 
are the generalization of the Witten-QCD model. They have found a butterfly shape and a double 
valuedness in the phase diagram of these confining geometries. In this paper, we additionally 
present the phase diagrams of the Maldacena-Nunez, Klebanov-Strassler and Klebanov-Tseytlin 
plus the Klebanov-Witten and Witten-QCD models by similarly calculating the length of the 
connected regions and the entanglement entropy of the connected and disconnect solutions. We 
also find the similar butterfly shape in the phase diagrams. We also find that if in a specific 
geometry the phase transition of Schwinger effect is fast and dramatic, this also would be the 
case in its entanglement entropy phase transition which is the case for WQCD and KT, and if the 
phase transition for Schwinger effect be mild, this also would be the case for its EE phase diagram 
which is the case for MN and KS. One can also compare the features with the conformal case of 
Klebanov-Witten and AdS as a limit of a mild transition. 

In section 6, we study the Schwinger effect in the presence of a magnetic field in addition to the 


electric field and study its effects on the pair creation rate. By adding a probe D8-brane, in 19 


and 1201 the authors studied the imaginary part of the Euler-Heisenberg effective Lagrangian, the 
rate of pair creation, the critical electric field and the effect of the parallel and perpendicular 
components of the magnetic field on the rate of pair creation in the background of Sakai-Sugimoto 
and the deformed Sakai-Sugimoto models. Similarly we calculate the DBI action and the imaginary 
part of the Euler-Heisenberg effective Lagrangian in our geometries. We find that in all of our 
supergravity confining geometries the parallel magnetic field would increase the rate of pair creation 
while the perpendicular magnetic field would decrease it which should be a general feature of all 
confining backgrounds. We then conclude with a discussion in section 7. 


2 The string profile in confining geometries 

In this section we present the PDE equations of the string profiles which can be used in finding 
the entanglement entropy of accelerating quarks moving on specific trajectories. 

As in |8||, by starting from the Nambu-Goto action for the AdSs geometry, 

ds 2 = R 2 [—^— u 2 dt 2 + u 2 (dx 2 + dy 2 + dz 2 )\, (2.1) 

u- 

and then by assuming the static gauge of (r, a) = (t , u ) and the embedding coordinate of X 11 = 
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(t,u,x(t,u), 0,0), one can find the determinant of the induced metric as 


\f—g = R 2 V 1 — x 2 + n 4 x' 2 , 


( 2 . 2 ) 


and the equation of motion as 


d (vAx' \ d f x \ 
du \y/^g) dt \^g) 


(2.3) 


Solving these PDE equations, in general is a difficult task, but for the case of AdS which enjoys a 
lot of symmetries and has a simpler form, the author in [8| could End the solution and therefore 
the string profile as a function of t and u as 

x=±\t 2 + b 2 —(2.4) 

V w 

If a heavy quark and an anti-quark accelerate on a specific trajectory due to the potential of an 
electric held, the classical solution from the Nambu-Goto action is a world sheet that is a part of 
AdS 2 and the locus of 


u 2 + (x 1 ) 2 — ( x °) 2 


M 2 
~E 2 ' 


(2.5) 


where M is the mass of the quark and E is the electric held. From this relation, the world sheet 
event horizon can be read as ue = yf- As in |4j, for a specihc trajectory such as a hyperbola, one 
can read the metric near the quark trajectory. 

Alternatively for Ending the induced metric on the world sheet, similar to |3j, one can use the 
Mikhailov relation between the embedding coordinate X m (t,u ) and the boundary quark position 
x m (t), (9], as 


X M (r, u) = ux^{t) + x m (t), X m (t,u) = (X^(r, u), u), (2.6) 

Then as in [5] one can End the proper area between the probe brane and the event horizon 
which is proportional to the free energy. By knowing the Unruh temperature that the quarks 
would feel in the accelerated reference frame in AdS as Tu = and by using the first law of 
thermodynamics one can read the entropy. If one can assume the semi-classical and heavy quark 
limit in the problem, all of those entanglement would be due to the “entanglement entropy” and 
so for the AdS the EE where found to be s = \/A, |4| where A is the ’t Hooft coupling constant. 

Now we look at the equations for the Af = 1 supergravity solutions that are dual to the 
confining geometries. One should notice that for this calculation, one first needs the string world 
sheet profile, the induced metric and the Unruh temperature of the accelerating particles in each 
background which we don’t present here. 
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2.1 Witten-QCD 


Between all the three mentioned confining geometries, the Witten QCD model is the most similar 


background to the AdS metric. In the string frame, its metric and dilaton field are 21 
ds 2 = (^) 3/2 (r] lw d/x ,L d/x v + ^-f(u)d9 2 \ + (—\ + R 3/2 u 1/2 dQl, 

R V 9^o / \UJ f(u) 


Ur 


m = i--^, 


u c 


R = (jtN g s )z a' 2 , 


e — 9t 


u 


3/4 


A 3 / 4 ' 


(2.7) 


If we consider the static gauge, and the embedding coordinate as = (£, w, #(£, w), 0, 0, 9{t, u ), 0, 0, 0, 0), 
then the equations would be complicated. So we assume 6 is a constant, and therefore the deter¬ 
minant of the induced metric simplifies to 





which leads to the following equation of motion, 

/m 3 — d / u 3 x'\ d f x \ 
V R 3 u 3 ) du \^g) dt \V~gJ 


( 2 . 8 ) 


(2.9) 


which is quit similar to the AdS case. If one finds the analytical solution of this PDE equation, 
similar to the AdS case, one can follow the procedures of |4j and find the entanglement entropy of 
heavy accelerating quarks in this model. 


2.2 Maldacena-Nunez 


The MN metric is obtained by a large number of D5 branes wrapping on S 2 , 22 . In the string 


frame the metric and the fields are 23 


ds 2 w = e^l—dt 2 + dx\ 2 + dx 2 + dx 3 2 + e 2h ^ r \d0i 2 + sin 2 didcfti 2 ) + dr 2 + ~(w l — A 1 ) 2 ], (2.10) 


where 


A 1 = — a(r)d9i, A 2 = a(r) sin 9id(j>i, A 3 = — cos9id<f>i, 


( 2 . 11 ) 


and the of ’s parameterize the compactification 3-sphere which are 


oj 1 = cos il>d9 2 + sin -0 sin 9 2 d(j) 2 , oj 2 = — sin r ipd9 2 + cos -0 sin 9 2 d(f) 2 

oj 3 = d'l/j + cos 9 2 d(j) 2 , 


( 2 . 12 ) 
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and also the other parameters of the metric are 


a(r) = 


2 r 


sinh 2 r 


e 2/ ' = r coth 2r — 


r 2 1 9 . 9 . 2e /l 

_ p -20 _ -20o_ 

sinh 2r 2 4 ’ 


sinh 2r 


(2.13) 


For the case of Maldacena-Nunez model, we assume the embedding coordinate as (t,r,x(t,r)) 
and all other coordinates will set to be zero. Then one would get 


a/-# = (1 - x 2 )(l + x /2 ). 

So the equation of motion of the string profile in this background is 

d fe 2<t> x'(l — x 2 )\ d ( e^xil + x' 2 ) 
dr 


(2.14) 


V~9 


dt l (x' 2 + 1)(7 2 — 1) 


= 0. 


(2.15) 


Again by solving this equation analytically one can find the string profile and then the entan¬ 
glement entropy of heavy accelerating quarks. 


2.3 Klebanov-Tseytlin 

The Klebanov-Tseytin metric is a singular solution which is dual to the chirally symmetric phase 


of the Klebanov-Strassler model which has D3 brane charges that dissolve in the flux 24 
The metric is 


ds 2 l0 = h(r ) ly/2 [—dt 2 + dx 2 ] + h(r) 1//2 [dr 2 + r 2 ds? 1 , 1 ] . (2.16) 

Here 7s^ M is a base of a cone with the definition of 

rf4M = i(3 5 ) 2 + ly(9‘) 2 . (2.17) 

i =1 

It is the metric on the coset space T 1,1 = (£77(2) x SU(2))/U{1). Also g l are some functions of 
the angles 0\, 9 2 , 0i, 02 ,0 as 

g 1 = (— sin 6id(j)i — cos 0 sin d 2 d(j) 2 + sin 077 2 ) / \/2, g 2 = (dd i — sin 0 sin 0 2 d(j) 2 — cos 0702) / x/2, 
g 3 = (—sin 0(701 + cos 0 sin 0 2 d(p 2 — sin077 2 ) g 4 = (7$i + sin 0 sin 0 2 d(f)2 + cos 077 2 )/\/2, 

g 5 = 70 + cos7i70i + cos d 2 d<p 2 , (2-18) 

and also 

r 4 „ oi 

MO = -4 In-’ L* = -g s M 2 e\ (2.19) 

II S " 
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In this frame, the asymptotic flat region has been eliminated. Also r = r s is where the naked 
singularity is located. 

For the case of KT the embedding would be (t, r, x(t, r)) with all other coordinates zero. Then 


and the equation of motion of the string prohle would be 

d_ _ l ( x(^ + fe") \ = 

dr V HV~9 ) dty ) 


( 2 . 20 ) 


( 2 . 21 ) 


2.4 Klebanov-Strassler 


The Klebanov-Strassler (KS) metric which is known also as warped deformed conifold is obtained 


by a collection of N regular and M fractional D3-branes 25 
The metric is 


ds \o = h 2 (^dx^dx^ + h 2 (r)dse 2 , 


( 2 . 22 ) 


and dsg is the metric of the deformed conifold which is 


ds 6 2 = —e 3 K (t ) 


LSi + (9 ' )2) + COsh2( i>l (93)2 + (9l)2 l + Binh2 (|M(9p 2 + (9 2 ) 2 ] 


(2.23) 


The parameters of the metric are 

(sinh(2r) — 2 r)s 


K(t) = 


h(r) = 


2 3 sinhr 

. . , x coth x — 1 , . . 

I(r) = / dx -g-(smh(2x) — 2x ) 3 , 

J T sinh x 


(2.24) 


and 


g 1 = — -j=[— sin6id(j)i — cos -0 sin d 2 db 2 + sin^d0 2 ], g 2 = —^=[d$i — sin -0 sin $ 2 d 02 — cos^d^]. 


v /2 


\/2 


g 3 = —=[— sindid 0 i + cos^ sin 6 l 2 d 02 — sin^d^], g 4 = —p[d$i + sin xjj sin d 2 d </>2 + cos^d 0 2 ], 


v/2 


a/2 

g 5 = d'l/j + cosdid^i + cos 0 2 d</> 2 . 


(2.25) 


For the case of Klebanov-Strassler similar to the KT, the embedding would be (t,r,x(t,r)) 
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with all other coordinate zero. Then 


= 


\ 


r /2 


(1 - i2) lT + « 4 ) 


(2.26) 


and the equation of motion of the string profile would be 


d fx'(l — x 2 ) 
dr V hyf=g 


d 


/2 

X ,X 


4 

63 


dt l \/—g h 6A' 2 (r)' 


= 0 . 


(2.27) 


2.5 Klebanov-Witten 

The Klebanov-Witten solution is similar to the KT throat solution but with no logarithmic warping 


26 . Unlike the other four mentioned metrics, it is a conformal non-confining geometry. We study 


this background to compare our results with the conformal case. 
The metric is 


ds 2 = h 2 g IJV dx ,J 'dx 1 ' + h' 1 ( dr 2 + r 2 ds^ 1 , 1 ), 


where, 

h = — t, and L 4 = g s N(a ') 2 . 

r 4 4 

The embedding would be (t,r,x(t,r)) with all other coordinates zero. Then 

\f zr g = ^/(l-i 2 )(l + ^-), 

and the equation of motion of the string prohle is 


d (x'(l-x 2 )\ d fx (l + f)\ 
dr V h^-g ) dt [ ) ~ " 


(2.28) 


(2.29) 


(2.30) 


(2.31) 


Therefore again solving this equation would give the string prohle and the induced metric near 
the accelerating particles’ trajectory in the background of the conformal KW model. 













3 The free energy of accelerating qq in the Miknowski 
background 

In [4] by minimizing the Nambu-Goto action and by using the solution of the PDE for accelerating 
particles in AdSs, the authors found the metric near the quark and anti-quark trajectory as 


ds 2 = \f\a! 


1 E 2 \ 2 2 ' 

V‘~Mi dT ~V‘ dTdu 


(3,1) 


then calculating the proper area between the probe brane and the event horizon yields 

Sn = —2 


M - eM - — 
2 


Tp. 


(3.2) 


Knowing the Unruh temperature of AdS, Tjj = and the free energy, = M — 4^ — y/XTu, 
the entanglement entropy of the accelerating quark and anti-quark have found to be s = y/X. 

Now the flat geometry can be the UV limit of the Hard-Wall and Witten-QCD model. Therefore 
for the flat metric 


ds 2 = y/~\a\—dt 2 + dx\ + dx\ + dx\ + dx\ + du 2 ), 
one can repeat the calculation of Semenoff and Hubeny. So one would have 

TVt — P 1. — 1 + x ' 2 , Itu — 


(3.3) 


(3.4) 


and then the equation of motion is 


d u 


x 


\/l — x 2 + X‘ 


a 


-d t 


x 


\/l — X 2 + X 


12 


= 0 , 


(3.5) 


and the solution of this PDE can be found as x(t,u ) = y/t 2 + b 2 — u 2 , where as in the previous 
case, for the accelerating quark and anti-quark, the constant is b = A-. From this solution one can 
see that the world sheet event horizon is at Ue = yy Now by using this solution, the components 
of the induced metric can be found as 

u 2 -b 2 _ t 2 + b 2 _ t 2 + b 2 

7tt “ t 2 + b 2 - u 2 ’ luu ~ t 2 + b 2 -u 2 ’ 1tu ~ t 2 + b 2 — u 2 ' ^ ’ 

Similar to the conditions that have been applied to derive the induced metric of Semenoff’ 
work, one can similarly reach to the following induced metric, 


ds 2 



2 E 2 

~w 


u 2 drdu 


(3.7) 


9 



















The D3-brane is located at % = Yrr. So the proper area between the probe brane and the 


event horizon is 




, , UE E 2 u 2 , M-, 

dr I - du + —J 


a j u „ M 2 

^2 JU M 


n r Vx f M E 2 y/X o \ M-, 

^ Tp [~ 67 ( ¥ “ M 2 27tM ' + ^‘ 


(3.8) 


Assuming that the Unruh temperature of the accelerated frame is Tjj 
would be, 


27rM 


VX t 2 a 2 


M 

2 V 12tt 2 T ' 12vr 2 M 3 + T 


then the free energy 


(3,9) 


So knowing the solution of any of the above PDE equations in any confining geometries can 
similarly lead to the proper area between the probe brane and the event horizon and yields the 
free energy of the accelerating qq in those geometries. In a similar way the holographic Schwinger 


effect were also studied in other geometries such as in de sitter space 27 . 


4 Potential analysis of the confining geometries 

One would think that for searching for any possible relationship between the entanglement entropy 
and the Schwinger pair creation rate, it would be interesting to first find the phase diagrams and 
the phase transitions for both quantities in few different confining backgrounds. 


In 10 the authors demonstrated some generals features in the phase diagram of Schwinger 


effect in all confining geometries and then in 11 , they have studied the “AdS soliton geometry 


as a special case. In these papers, the authors studied the potential of a confining background and 
then they have plotted the total potential, Vtot(^) versus the distance between the quark and anti 
quark x. Their setup and the world sheet configuration for the quark and anti quark potential is 
also shown, in Figures [l] and [2j To compare with our diagrams, their plot of the phase diagram 
for the AdS soliton geometry [l0| is re-produced in Figure |4j 

Now in the following sections, we do a similar calculation for our class of supergravity confining 
geometries and then numerically we find the phase diagrams. 


4.1 Witten-QCD 

First for the Witten-QCD model, if we assume that a probe D3-brane is located at u = u 0 and 
then we assume the following ansatz for the Wilson loop, 


10 











probe 



x 

A 


r t r c 


h) 


Boundary 


► r 


Figure 1: World sheet configuration in 3D 


Figure 2: World sheet configuration in 2D 


x° = r, x 1 = a, u = u(a), 6 = 9(a), 


(4.1) 


the NG action would be 


C = v /= ^ = 




u' 2 fR 


1 rx \U 

1 7,3 ' 


(4.2) 


As the Lagrangian does not depend on a , the following Hamiltonians are conserved, 

u _ dC x r u _ dC a a r 

o(d a u ) o(d a 6) 

Then there should exist two different minimal surfaces that satisfy the following relations, 


(4.3) 


5^ = 0, at u = u c ( u t <u c < u 0 ), ‘ ! f = I), at 0 = 9 C ( O t < O c < 9 0 ). (4.4) 

aa aa 

However by assuming 6' = 0, we consider a rectangular Wilson loop which makes the calculation 
simpler. So 


£ = 


\ 


1 / du\ u 3 

1-4 \ dcr ) + R3 ’ 


(4.5) 
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and then from the conservation of H u we would get, 


du 

da 



(4.6) 


By integrating the above equation one could get the length of the string between the quark and 
anti quark in WQCD as 


x = 2 i?a 




Now by defining the following dimensionless quantities, 


one can simplify x as 


u 

u c 


•) 

a, — —, 

b = — 

He 

U 0 

u 0 


2 Rl 


x = 


dy 


(uqO) 


2 J 1 


y fe 3_l)(y3_(|)3 


Then the sum of the potential and static energy would be 


If 


Vpe+se — 2 Tp / daC — 2TfU 0 a / dy — -j= - 

J ° h -!)(>-£ 


(4.7) 


(4.8) 


(4.9) 


(4.10) 


For the large x limit, (a —> b), the sum of the potential and static energy is 

3 

Vpe+se = Tp - ^ 3 ^ x + 2 TpUob — 1 ^ . ( 4 - 11 ) 

So from the first term which is the quark and anti-quark potential, we can read the confining string 
tension as 

°st= TF Qy. ( 4 . 12 ) 

This matches with the result coming from the relation a s t = • Also the second term gives the 

static mass of the quark and anti-quark, 


2 Tp(uo — u t ) = 2 mw- 


(4-13) 
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Then from the DBI action, one can read the critical electric field for this geometry as 


E c 


MS' 


(4.14) 


Now we can define the dimensionless parameter a = So the total potential energy would be 

-C'c 


2T F u 0 a 


Vtot = Vpe+SE — Ex = 


1 

a 


dy 


vV- i)fe 3 -5) 


2T f uqC( 

I 

a 2 



_ dy _ 

yV - i)(y 3 - §) 


(4.15) 


Vx_tot 



Figure 3: The plot of total potential versus x for Witten QCD model when b = 0.5 and 2 Tp = uq = 1 . 


By comparing the diagrams below, one can see that the form of plots are similar with slight 
differences in small x limit. One can see that in both of them, there exist three phases, one with 
stable potentials and no pair creation, one with exponentially suppressed potentials with tunneling 
pair creation and one with catastrophically unstable potentials with exponential pair creation. 


As it has been demonstrated analytically for a general background in 10 this behavior is 
universal in all confining geometries. However there are still some minor differences between the 
phase diagrams of different confining backgrounds which here we aim to detect and then compare 
with the phase diagrams of the entanglement entropy of these backgrounds. 

From these plots one can see that in Witten QCD diagram, for a = 0.01 there is no zero other 
than the origin and no Schwinger effect can occur. For a = 0.1 the potential becomes flat. For 
a = 0.25 or 0.6, there is a barrier in the potential which as it can be seen is different from the AdS 
soliton case with a bigger curvature. In this phase the Schwinger pair creation can only occur by 
tunneling through this barrier and the rate of pair creation is strongly suppressed. For larger E 
as a = 1.3 for instance, the potential becomes catastrophically unstable and the Schwinger effect 
would occur and therefore a current can be created. In this phase the probability of pair creation 
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Vx_tot 


Vx_tot 




Figure 4: AdS soliton Figure 5: Witten QCD 

would not be any more exponentially suppressed. The curvatures are bigger for the Witten QCD 
case and a small bump can be seen at around x = 0.5 which is not present for the AdS soliton case. 
Also comparing all the diagrams one can see that in Witten-QCD and KT the phase transition 
happens faster and in a more dramatic way relative to the other backgrounds. 


4.2 Maldacena-Nunez 

Now we repeat this calculation for the Maldacena-Nunez background. We assume x° = r, x 1 = 
cr, r = r(a), so the determinant of the induced metric and therefore the Lagrangian is £ = 
e4>yJ r i2 _|_ i Since 

dC 

———d ff r-£, (4.16) 

o(d a r ) 

is a constant, we would get 


r 


/ 



(4.17) 


So, £ = 4-. Also from the DBI action one can find the critical electric held as 

' p<Pc 


e<Pc 


E r = Tl 


g2 4>{ro) 
e <K r c) 


(4.18) 


We assume that the probe brane is located at r — r c — 1, and for the sake of similarity to 
the previous calculations we define a = W therefore r 0 = A One should also notice that the 
MN geometry ends when becomes undefined which happens at the roots of e h = r coth(2r) — 
sinh[2r 2 ] — 4 w ^ iere plot is shown in Figure J^J. 

One can see that the positive root is approximately located at r t = 0.6 where the geometry 
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Figure 6: The Plot of e h versus r . 


ends. Now by defining 


x = 



dr 



J 


and the total potential as, 


Kat 




(4.19) 


(4.20) 


one can find the phase diagram in Figure [7} 



Figure 7: The plot of total potential versus x for the Maldacena-Nunez model. Here do = 0, r'o = Tp = 1. 

By expanding the potential Vpe+se for large x, i.e, r —>• r c , one can find the string tension and 
quark mass as, 


a st = e^Tp, 2 m w = 2Tpe^ c (r 0 - r t ). 


(4.21) 


This result for the string tension matches with other relation from One can see that still 

there exists three different phases. For a till around 0.47 no Schwinger effect would happen. For a 


















between 0.47 and 0.48 Schwinger effect would occur as a tunneling process, and for a larger than 
0.48, the potential becomes unstable. 


4.3 Klebanov-Strassler 


Now for the KS background we assume x° = t, x 1 = a, t = r(c r). So the components of the 
induced metric are 


7 tt = -h 2(r), 


lo 


= h 


r) + 


T 


/2 /r5 (r)e 


6/1 2 (r) 


(4.22) 


Therefore 


C = \ h~ l {r) + 


/2 A 
t z e 3 


6/1 2 (r) 


Again knowing that — £ is constant, one can derive r' as 


/ \/6 A (r) fT7 — - 

T = — TttVW - KV- 
6 3 h(r) 


Also the critical electric held is 


E x = T F h 2 (tq). 


(4.23) 


(4.24) 


(4.25) 


At r = r c one would have t' c = 0. We assume To = 1, so one can find 

1 


x = 2 


63 h(ya) 


V6K(ya ) ^(a) - h(|/a) ’ 


and the total potential as 


‘lETp r a yjh(a) 1 

i/6 Ji K(ya) yjh{a) — h(ya ) 

2 eiTpoi h(ya) 1) 

\/6 J \ K (ya) yjh(a) — h(ya,) 


(4.26) 


(4.27) 


The plot of KS phases for different ct is shown in Figure [8} 

The numerical calculation of this case would take more time to finish but the final phase 
diagram in general is similar to the three other confining geometries. The part that are missing in 
the plot are due to some singularities in the numerical calculations of the integrals. Again three 
phases for different a can be detected and the phase transitions are mild. 
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Vtot 



Figure 8: The plot of total potential versus x for the Klebanov-Strassler model, assuming tq = (3 = 
Tp = e = 1 . 

4.4 Klebanov-Tseytlin 

The calculation can be repeated for the KT background. Taking x° = t, x 1 = er, r = r(cr), the 
components of the induced metric are 7 t t = h~^, = —(h~^ + r' 2 h 5 ) and so £ = \/h~ l + r' 2 . 

Again using the conservation of the Hamiltonian and knowing that at r = r c , r' = 0, then we 
find 


./ = y /h ( r c) ~ h(r) 
h(r) 


(4.28) 


After taking the integral of this relation and by defining y = (j 1 , a = ^ and b — one can 
find 


x = 9v^2Me 2 y/~g~ s 


dy 


1 »V hl (f)-i? hl (¥)’ 

Vpe+se = 2T f [ Cda — 2Tparn, \n - f — - 

V bJi , /in | i In fi 


(4.29) 


From the DBI action the critical electric field is 




L 2 


(4.30) 
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By defining a — Jr-, the total potential is, 


Vtot = 2T F ar 0A /ln- J 


dy 


2aT F rl 


dy 


In ? - 4 In 7- 


y 4 b 


Ini J ‘ 


(4.31) 


Now by simplifying this relation by assuming b = 0.4 and ro = M = T F = 1, one can find the 
plot of Vtot versus x which is shown in Figure [9j 



Figure 9: The plot of total potential versus x for the Klebanov-Tseytlin model for b = 0.5 and ro = 
M = T f = 1 . 


The Klebanov-Tseytlin model is a limit of Klebanov-Strassler and there is no wonder that the 
phase diagram would be very similar. Again three phases can be seen for different a. However in 
these two cases the exact numerical values of a does not correspond to each other, since we didn’t 
match the numerical constants. 


4.5 Klebanov-Witten 

Now for comparing our results of the phase diagrams of confining geometries with the conformal 
backgrounds, we also study the Klebanov-Witten geometry. The phase diagram would be exactly 
identical to the phase diagram of AdS case which has been studied in 12 . Assuming x° — t, x 1 = 
a, r = r(a), where h — Then £ = \fh~ x + r' 2 , and as the Hamiltonian is conserved, one would 
get, 


r 


/ 



(4.32) 


where ^ = 0 at r = r c . So £ = and from the DBI action the critical electric field is E x = -jy? 2 -. 
By integrating r' and by defining y = a = ^ and b — (r t — 0 here) the distance between 
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the quark and anti-quark and the potential are 


2 L 2 [<> dy 

T — / 

r o a J i y 2 \Jy i — l' 

(4.33) 

Vpe+se - 2 T F r 0 a j - 

Ji vV - 1 

(4.34) 


So for large x, as a —> b = 0, the potential is “zero” which is a feature of non-confining, conformal 
backgrounds. Also, 2 mw = 2T F r 0 (l — a) and for a = b = 0, it gives 2 mw = 2 T F r 0 . Then 


Vtot = 2 T F r 0 a 



V 2 dy 

V 7 ?/ 4 -! 


2 T F r 0 a 
a 


dy 

y 2 Vy 4 ~ 1 


The phase diagram of KW and AdS are exactly similar and is shown in Figure 10 


V tot 



(4.35) 


Unlike the other four confining geometries, as the KW and AdS are conformal geometries, there 
are only two phases here. Even for a very small a or electric fields, there is a zero at larger x and 
so the pair creation would happen by a tunneling process. The other phase is the unstable one for 
a larger a. 


5 Entanglement entropy of a strip in confining geometries 

Now we would like to study the phase diagrams of entanglement entropy in these confining ge¬ 
ometries and then compare the phase transition in EE with the Schwinger effect and compare 
the diagrams of different geometries. So one way is to use the method in |4| for calculating the 
entanglement entropy of accelerating quark and anti-quark. But as we have mentioned in section 
2, one first needs to solve several PDE equations analytically. Other methods could be using the 
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ideas in 281 or 29 , to study the entanglement entropy of local operators or localizes excited states 


in each background. But here we follow the calculation of 17 and 16 to find the entanglement 


entropy of a strip in each confining geometries. Also the entanglement entropy of multiple strips 


can be calculated similar to 30 


So as in 16 , based on KKM, the generalization of Ryu-Takayanagi conjecture for the non- 


conformal theories would be 


S = 


G 


( 10 ) 

N 


d 8 ae 


G 


( 8 ) 

ind’ 


(5.1) 


The authors showed that the entanglement entropy can be found by minimizing this action over all 
surfaces that ends on the boundary of the entangling surface and there are actually two solutions, 
where one of them corresponds to a disconnected region and the other one is a connected surface. 


In 29 , the authors showed that always, only one of the two possible configuration would dominate. 


So if one writes the gravitational background in the following form of, 

ds 2 = a(p)[/3(p)dp 2 + dx^dx^l + gijdd l dd^, (p — 0,1, ...d), (i — d — 2,..., 9), (5.2) 


then the volume of the internal manifold is Vjnt = / d0 A/set[(pj], and one can define another 
function H(p) as 


H(p) = 


a 


(5.3) 


This function for confining geometries is a monotonically increasing function while /3(p) is a mono- 


tonically decreasing function. Using the KKM equation 1161, 17 the EE for the connected solution 


is 


o,.' Vd-i r, i mh(p) 

(io) / dP< . _ H(n o) 
j po y i tf( p ) 


(5.4) 


and the EE of the disconnected solution is given by 


Vd i 


SdM = / dp^MH(p), 

2 G n J Pa 


(5.5) 


and the length of the line segment of the connected solution is 


L(Po) = 2l dpl-g^ 

p0 V H(po) 


- 1 


(5.6) 
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The difference of the connected and disconnected solution would be finite, so S is defined as 


S(po) = 


2 G 


( 10 ) 

N 


Ki-i 


(S c — S D ) = / dp. 


' P 0 


I MH(p) 

i _ H(po) 
H{ P ) 


dp\ZP(p)H{p). 


(5.7) 


'PA 


Now we study L and S for our specific geometries. 


5.1 Witten-QCD background 


For the case of Witten QCD background similar to 
and H{u) as 


17 , one can define the functions a(u),/3(u) 


a(u) = 




( 87 t 2 \ 2 AR 6 u 5 f(u) 

\ 3 / 9u 0 gi 


(5.8) 


where f{u) = 1 — From the functions of j3 and H, it can be seen that as it has been suspected, 
P(u) is monotonically decreasing, H is monotonically increasing, H shrinks to zero at u = 0 and 
/3(u) diverges at u — u t — 1. 

The EE for the connected and disconnected solution are respectively 


Sc(u 0 ) = 


E 2 8n 2 R2 
G n ^ 9 g 2 sy /u t , 


du 


UQ 


\ 


u z 


i _ u o f(u o) ’ 
u 5 f(u ) 


Sd(uo) — 


V 2 87 r 2 Ri 


udu. 


G^9g 2 sy /u t J UA 

Also the length of the line segment of the connected solution is 


(5.9) 


L(u 0 ) = 2 


du 


'UQ 


\ 


(R)3 _J_ 

W ' /(«) 


f(u)u 5 

/(«o)«o 


- 1 


(5.10) 


Their plots are shown in Figures [111 and [l2| 

Here we take u t — u\ — 1. Due to the peak in the plot of L and also the butterfly shape 
and the double valuedness in the plot of S, one can deduce that a phase transition and therefore 
a confinement phase exist. This phase transition and the shape of the peak are similar to the 
instability of V and the phase transition in the previous section. One can see that the phase 
transition is more dramatic relative to the other geometries. This dramatic phase transition 
in smaller x were also seen in the Schwinger phase transition in WQCD relative to the other 
geometries. So there might be a deeper relationship between these two quantities. 
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Figure 11: Plot of L(uq) v.s uq for WQCD model. 


Figure 12: Plot of S vs. uo for WQCD model. The blue 
line is the connected solution and the dashed red line is 
the disconnected solution. 


5.2 Klebanov-Tseytlin 

For the KT case, the volume of the internal part, ds T i,i, is Eyb So for the KT background the 
functions would be 

1 1 f)7T^ ^ f 16 \ 2 

a(r) = h(r)~ 2 , /3(r) = h(r), V) n t = -^~r 5 h(r)*, H(r) = ( — ) 7r 6 r 10 /r(r). (5.11) 

Again f3 and H show the expected monotonically decreasing and increasing behaviors. Then 
the length of the connected region and the EE are 


L(r 0 ) = 9V2My/g~ s e 2 / dr 


'r 0 


Sc(ro) = 


r»(io) / 

U N Jr 0 



lr* 1 

ro 6 In — 1 

V 

r s 

) 

r In — 

dr- 

r s 


, 6 In 31 
Y _ Lo _is 


r 6 In , 

r 


c , , 121 WM*g,e‘ f' , , . 

SdM = - - / *’■>%- 


G 


(10) 

JV 


(5.12) 


'r A 


The plots for the M = e = g s = 1 and r s = 0.5 is shown bellow in Figures 13 and 14 


As one can see the form of the plots are very similar to the KS geometry if the naked singularity 
is placed at small r as in here it is set to r s = 0.5. However increasing r s causes that the red dashed 
line goes down and the butterfly shape of the diagram get a flatter curvature as it is shown in 
Figures [l3j [l4j [15] and [16 
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Figure 13: Plot of L(ro) vs ro for KT, r s = 0.5. 
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Figure 14: Plot of S(ro) vs L for KT, r s = 0.5. 
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Figure 15: Plot of L(r 0 ) vs. r 0 for KT, r s = 2. Figure 16: Plot of 5(r 0 ) vs. L(r 0 ) for KT, r s = 2. 


5.3 Klebanov-Strassler 

For the KS background, the functions are 


a = h a(r), P(t) = 


8vr 6 


h{r)e 


4/3 


6 K 2 (t) ’ 


H(t) = — e 20/6 h(r)A' 2 (r)sinh 4 (r). 
o 


(5.13) 
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Again, (3 is a monotonically decreasing function and H is a monotonically increasing function. 
Now using f3 and H we study the entanglement entropy of a strip in this geometry. So 


L(r 0 ) = 


2 let 

"W 


dr- 


sinh(r) 


'TO 


(sinh( 2 r) — 2 r) 3 


\ 


h(r) 


sinh(r) 2 f sinh(2r)—2r ^ 3 
sinh(ro) 2 V sinh(2ro)— 2tq 


Sc(t 0 ) = 


V 2 tt e 


3 A 


3 G 


( 10 ) 

N 


dr 


h(r) sinh 2 (r) 


'TO 


sinh(rp) \ 
sinh(r) I 


sinh(2ro)—2 tq \ 3 


Sd(t, q) — 




3 ,4 


3G 


( 10 ) 

N 


sinh(2r)—2r 

drh(r) sinh 4 r. 


f T A 


(5.14) 


L(r 0 ) 



Figure 17: Plot of L(uq) vs. u for KS. 


S(r 0 ) 



Figure 18: Plot of S{uq) vs.uq for KS. 


The plot of the length of the connected solution, L and the entanglement entropy, S is shown in 
Figures [TT] and [18} The entanglement entropy of KS case with dynamical flavors were also studied 
in |31|. From the behavior of L and the butterfly shape of S one can detect the confining phase. 
One can see however that the phase transition is milder relative to the Witten QCD and KT cases. 
This milder phase transition were also seen for the Schwinger effect phase transition. 


5.4 Maldacena-Nunez background 

For the Maldacena-Nunez background, the functions are 

a(r) = e^, /3(f) = 1, V int = 8n 3 e 2h e^, H(r) — An 6 e 8 ^ 0 (sinh. 2r) 4 . (5.15) 

Unlike the other backgrounds that we study here, for the Maldacena-Nunez metric, the function 
f3 is a constant and is not monotonically decreasing. Also H(r ) is not monotonically increasing. 
The L(tq) and EE functionals are 


24 




















2 


Sc(r 0 ) = 


L(r 0 ) = 

V 2 TT 3 e^° 


dr 


'r o 


(sinh 2r) 4 


- 1 


G 


( 10 ) 


r o 


S D (r 0 ) = 


N 

V 2 i r 3 e 4 ^ 0 


G 


( 10 ) 

N 


(sinh 2ro) 4 

00 , (sinh2r) 2 
dv — — 

-i (sinh2ro) 4 
1 (sinh 2r) 4 

/ dr (sinh 2r) 2 . 

'r 0 


Their plots are shown in Figures |19] and 20 


(5.16) 
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Figure 19: Plot of L(ro) v.s ro for MN model. 


r ° * 0.7848 ' 0.7850 0.7852 0.7854^^°^ 

Figure 20: Plot of S vs. L(ro) for MN model. 


5.5 Klebanov-Witten 


As the Klebanov-Witten geometry is not confining, it would be interesting to compare the behavior 
of the functions f5 and H and also L and S with the confining geometries studied above. 

As one can see in figures |2l| and 22 there is no phase transition in the plot of S and the true 
solution is the connected one. Also there is no peak in the plot of L which again specifies that 


KW is indeed a conformal geometry. This can be compared with the diagram of 10 
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T7 ,. 01 loi e t c r/117 Figure 22 : Plot of S vs. L for KW. 

Figure 21 : Plot of L vs. p for KW. 


6 The critical electric field in the presence of magnetic 
field 

Now in this section we assume that in addition to the electric held, a parallel and a perpendicular 
magnetic held components are also present. By using the Euler-Heisenberg Lagrangian, we then 
study the critical electric held which would lead to the Schwinger pair creation in our four confining 
geometries. We see that similar to the Sakai-Sugimoto and deformed Sakai-Sugimoto models, the 
parallel component would increase the pair creation rate and the perpendicular component would 
decrease it. 


6.1 Maldacena-Nunez 

For the MN metric, we have 


?<t> 0 


^string—MN 


yj sinh (2r 


27ra'v / 2 /V*coth (2r* 


sinh-(2r*) 




( 6 . 1 ) 


Based on the assumption that a geometry has an IR wall and the probe D-branes would hit the wall, 
by calculating the imaginary part of the DBI action with a constant held strength, the authors 
in 20 showed that the critical electric held, where the Euler-Heisenberg Lagrangian becomes 


imaginary and therefore the Schwinger effect start to take place, is at E c which is, 


E r , = o, 


string \ 


rr 2 . I I 5|2 

u string T pi 

\ . 15 12’ 

^ u string ' 


( 6 , 2 ) 
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Here airing = 2 ^ anc ^ r * where the DBI action vanishes. By calculating the equation 
can get 


6.2 


one 


E c -mn — 


e^° y /sinh (2r* 


r>< 


27ra' (4 coth (2r*) r* — 1 — 4csch 2 (2r*) r 2 ) 4 
e 2 0 ° ( 2 r*) + 47 r 2 a , 2 H 2 4 coth (2r*) r* — 1 — 4csch 2 (2r*) 

^ g2(/> 0 si n h (2r*) + 47r 2 a /2 i?|| 2 ^4 coth (2r*) r* — 1 — 4csch 2 (2r*) r 2 ^ 


(6.3) 


For Hi 1 = H = 0 the critical E is 


Ecr-MN ~ 


e^° >/sinh ( 2 r*) 


(6.4) 


‘1'Ka! (4 coth (2r*) r* — 1 — 4csch 2 (2r*) r 2 ) 4 
We can re-derive this result by finding the imaginary part of the Euler-Heisenberg action 


and therefore we can find an indicator for the universality of the equation | 6 . 2 | for the confining 
geometries. 

We then calculate the DBI action by finding the induced metric of MN on D3-brane or D7-brane 
profile. The D3-brane profile is 


ds 2 = (1 + ~t) ^ (— dt 2 + dx 2 ) + (1 + ~~r ) ^ (dr 2 + r 2 dQ 5 2 ), 


(6.5) 


and the D7-brane profile is 


9 r 2 R 2 

ds w = -mdxpdx 1 * + —ds 2 e j, 


R2 ■ r z 

ds 2 6) = dr 2 + r —{^.(w l 2 + w 2 2 ) + \w 3 2 + (dd - ]-f 2 ) 2 + (sin 9d(j) - \fi) 2 ), 


3 v 4 v-, ■ — ' 1 3-- 1 2 J " / ' v r 2 

?re R 4 = ^TcggNcis and R is the AdS 5 radius. 

Now we calculate the DBI action by finding the induced metric on D7-brane. So 


( 6 . 6 ) 


ZZmn—D 7 — —T 7 I d 2 x dfl 3 E 0<> I dr 


00 coth( 2 r) 


r z _ 1 

sinh(2r) 4 


'Hck 


yj\ sinh( 2 r) 


■x 


2 e 4< N sinh 4 ( 2 r) 


(2vra) 2 e 2 ^ (F 0 i 2 - F 12 2 - F 13 2 - F 23 2 ) sinh(2r) 


^ 1 _ gr 2 - cosh(4r) + 4r sinh(4r) ^4r coth(2r) - 4r 2 csch 2 (2r) - 1 


(2t r «) 4 F 01 2 F 23 2 . 
(6.7) 
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For Foi = E ], F23 = Bi, F 13 = Bo, F\ 2 = B 3 we derive the Lagrangian as 


EmN-D 7 — —T7 I (Fx dh^i? ^ () 


r-oo (r coth(2r) 
dr —— 


sinh(2r) 


l 

1 \ 4 


'Hck 


sj\ sinh(2r) 


■x 


\ 


2 e 4< ^° sinh 4 (2r) 


(27ro;) 2 e 2<?;>0 (B 2 — Ei 2 ) sinh(2r) 


1 — 8 r 2 — cosh(4r) + 4 r sinh(4r) 


4r coth(2r) — 4r 2 csch 2 (2r) — 1 


- (2vra) 4 J B 1 2 J E 1 2 . (6.8) 


By solving the E which makes this Lagrangian imaginary, we find the same E as equation 6.2 So 
this way the universality of this relation can be checked again. From Figure [23] one can see that 


\1N Critical E 



Figure 23: The E cr versus parallel and perpendicular magnetic fields for a = 1, <^o = 0 and r = 2. 


the behavior of the critical E vurses the magnetic fields is very similar to the general figure as 
in 


20 


One can specifically check that for this background, by increasing the parallel magnetic field 
the imaginary part of the EH Lagrangian would increase which leads to an increase in the pair 
creation, but increasing the perpendicular magnetic field decreases the rate. 

Now we look at the electric field dependence. For the case of Hu = 0, we have 


ImC 


perp.B 


f 

' 0.05 


l 

1\ 3 


dr ( r coth(2r) — r 2 csch(2r) — - ) x 


4sinh 3 (2r) 


2(2vra) 2 (B w 2 - E\\ 2 ) 


\ 8r2 + cosll < 4r > - 4rsillh ( 4r > " 1 4r c°t h (2r) — 4r 2 csch 2 (2r) - 1 


(6.9) 


and for the case of B± = 0, we have 





















Figure 24: The ImL for the MN background, vs. parallel and perpendicular magnetic fields for a = 1, 
cj )o = 0 and E = 10. (We have normalized the fields.) 



Figure 25: The ItuCmn versus electric field and perpendicular magnetic fields for B || = 0, a = 1, (f>o = 0 
and r = 2 . 
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dr ( r coth(2r) — r 2 csch(2r) 
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2(2na) 4: B\\ 2 E\\‘ 

sinh(2r) 


4 e 4 </>o sinh 3 (2r) 


8 r 2 + cosh(4r) — 4r sinh(4r) — 1 


2 (2vra) 2 e 2 ^ [B\\ 2 - E\\ 2 ) 
4r coth(2r) — 4r 2 csch 2 (2r) 


( 6 . 10 ) 
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Figure 26: The IttiCmn versus electric field and parallel magnetic fields for B± = 0, a = 1, = 0 and 

r = 2. 


6.2 D3 probe brane in MN background 

The pullback of the Maldacena Nunez metric (which consists of N D5 brane) on D3-brane world- 
volume is e^rj tw . So the Lagrangian would be 

£mn-D3 = — ^3 J d 4 x 1 /(e _4 ^°(—(27ra) 4 e 4 ^ 0 F 0 i 2 F 2 3 2 + (—1 + 4r coth(2r))csch 2 (2r) — 4r 2 csch 4 (2r) — 
(27ra) 2 e 2<?i0 (F 01 2 — F 12 2 — F 13 2 — F 23 2 )csch(2r) y —1 + 4r coth(2r) — 4r 2 csch 2 (2r))). (6.11) 

There is no integral of “r” for the D3-brane case, so we can actually study the Schwinger effect in 
different r by inserting D3-brane in the geometry as the probe. 

First we study the potential in MN geometry when all the electric and magnetic fields are off 
which is shown in Figure [27} 



Figure 27: The potential on the D3 brane when all the fields are off. 

As one can see, all the D3 branes are being pulled to the tip of r = 0 when they are inserted 
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in MN and KS, due to the effect of F 5 and gravity. 

Now we investigate the potential when the fields are on. If we take F 0 1 = E \, F 12 = B s , F 2 3 = 
B\, F 13 = B 2 then the potential is 

1 / = {e- A<t>0 {-(2'KaYB l 2 e^°E l 2 + csch(2r)((2vra) 2 e 2 ^(F 1 2 + B 2 2 + B 3 2 - E x 2 )x 
\J — 1 + 4r(coth(2r) — rcsch 2 (2r)) + csch(2r)(—1 + 4r(coth(2r) — rcsch 2 (2r))))) j 2 . (6.12) 



Figure 28: The potential on the D3 brane for Eu = 10, = 0, a = 1, </>o = 0. 

This is shown in Figure [28] One can see that the potential is looking like two walls which have 
higher slopes near the IR region (r —y 0). In the UV (r —y 00 ) the potential is zero and has zero 
slope. Increasing E\\ would increase the slope of the potential and the pair creation rate. 



Figure 29: The imaginary part of Lagrangian related to the pair creation for E || = 10,R|| = 0,a = 

1 ) <t>o = 0 . 
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Again from Figure [29j one can notice that the pair creation happens with higher rate near the 
origin, at r —y 0, and at bigger r it would decrease. Increasing the perpendicular magnetic field 
would decrease the pair creation till make it zero at a specific value of the perpendicular magnetic 
hied at any r. 

Now we turn off the perpendicular magnetic held and study the effect of the parallel component 
of the magnetic held. As it can be seen from the following three hgures, parallel magnetic held 
generally increases the whole pair creation rate, however at B\\ = 0 the pair creation is mainly 
happening at small r. Increasing B\\ increases the pair creation in the UV and decreases the area 
which the pair creation is happening in the IR till makes it zero there, but in total, increasing B\\ 
would increase the imaginary part of Lagrangian and therefore the rate of the Schwinger effect. 



15 


Figure 30: The imaginary part of the Lagrangian for Eu = 10, B || = 0, a = 1, <f>Q = 0. 
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Figure 31: The imaginary part of Lagrangian for Eu = 10, B\\ = 3, a = 1, (j>o = 0. 

It would be interesting to notice that when Eu = 0 even with a strong Bu, no pair creation 
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would happen. Only at small r, in the IR region, increasing the perpendicular magnetic field would 
increase the potential. 

One should notice that heavier states with higher charges lie in the IR as the hadrons’ wave 
functions fall as r -A and A ~ J. This can be one reason that we see such a behavior in the IR. 
One should also note that Wilson loops also find it more favorable to lie at the end of the space 
in the IR. 

The plot of — Det(FMN + gMN ) vs. r and B±, for a constant E = 10, and zero parallel magnetic 


field, Figure 32 , shows that there exist a critical r, ( for Bj_ = 0 is around r = 1.2 ) where a hole 
is forming. For confining theories there is a critical r 0 where ^r| r _ r = 0 which for our specific 


geometry gives r = 1.118 which is close to what we have seen from the figures. Increasing E would 
increase the radius of this hole. 



6.3 Klebanov-Strassler 

The Klebanov-Strassler (KS) metric is 

ds \o = h~^{i-)dxij,dx^ + h^[r)dsQ 2 . (6.13) 


The form of this metric is also like 6.18 
string tension is 


so the D-branes hit the IR wall. The Klebanov-Strassler 


(^string—KS 


h z 1 (t») 
27ra' 


(6.14) 
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Using this and the equation 6.2, we can find the critical E as 




^ M r ) _ 

27TQ;' \ 1 + (27ra'i?|j) 2 h(r) 


(6.15) 


For B — Bn — 0, as it was obvious from E cr = f^4, we would have E cr = h 2 (r) 


1 

IS 

2 TTOi' 


The behavior of the function h(r) is shown in Figure 33 It shows that after r > 10 it is 


practically zero. So for showing the numerical plot we don’t actually need to take the integral of 
the Lagrangian to r = oo, but rather we take as r = 10. So the Lagrangian for the constant field 


/l(T) 



Figure 33: The behavior of the function h(r) vs. r for the Klebanov-Stressler metric, 
strength is 


£ks = —TV / d A x dh^e 



1 + (2tto) 2 (£|| 2 + B ± 2 - T|| 2 ) h{r) - (2iTa) 4 B pL 2 E ]] 2 h(Ty 


h{r) 


(6.16) 


Again as in Figure [34] one can check that for the KS background by increasing the parallel 
magnetic field the imaginary part of the Lagrangian would increase leading to an increase in pair 
creation and vice versa for the perpendicular magnetic field, so by increasing the perpendicular 
magnetic field, the imaginary part of the Lagrangian would decrease leading to a more stable phase 
with lower rate of pair creation. 


6.4 Witten-QCD 

The Witten-QCD metric is 




-f{u)d0 2 + R*u*<m\. 


(6.17) 


34 



















Figure 34: The Im Cks vs. B± and i?| of the KS background for r = 9.5, h(r) = 0.000798174, and 
£j| = 10, a = 1. 



Figure 35: The Im Cks vs. B± and E of the KS background for a = 1, r = 9.5, h(r) = 0.000798174. 

Since the generic form of the confining geometries is 

ds 2 = g(r)rj [lv dx^dx 1 ' + f{r)dr 2 + h(r) [internal space], (6.18) 

in the Witten QCD model, the internal geometry mix with the radial coordinate u, but still it is 
confining and the probe D-branes hit the IR wall. 

For WQCD, a st = 2777 which is consistent with our calculation of the potential in the 
previous sections. 
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Assuming Tp — ^7 from the above equation one would get 


( 6 . 

and when the magnetic held is off, this leads to the familiar result E c = cr st . Now the DBI action 
in D5-brane background including a constant electromagnetic held is 

Sgf = ~T 5 J d A xdue~^-det(P[g} ab + 2vr a'F ab ) = 

2 T 5 g s lu 6 + 4 (Bl + Bf - E 2 )n 2 R 3 u 3 a 2 - 16£j£% 4 f2 6 a 4 

3 Rz JuJl fK 3 u 0 

Again one can see that in this model too, by increasing the parallel magnetic held the imaginary 
part of the Lagrangian and therefore the rate of pair creation would increase, while increasing the 
perpendicular magnetic held would decrease it. 


( 6 . 20 ) 



E r = Tp l 

il 


x 3 

U* \ 2 


T|(f ) 3 + S 2 


\ 7?(t ) 3 + B\ 


2 1 


6.5 Klebanov-Tseytlin 

The critical electric held in the presence of a magnetic held in KT background would be 


E c = T f — J hi- 


T 2 F £ln^ + \B\i 


%ki r ^ + \B\ 


\ n L i In 


( 6 . 21 ) 


and the DBI action in D5-brane background is 


DBJ t 5 f°° /r 8 + 4(5 2 +5 2 -^ 2 )L 4 vr 2 r 4 a 2 ln^-16BjE 2 L 8 7r 4 a 4 (ln(^)) 2 

<6 ' 22) 


which again one can check that increasing the parallel magnetic held would increase the rate of 
pair creation and increasing the perpendicular magnetic held decreases the rate. 


7 Discussion 

In this work we studied the Schwinger effect and the entanglement entropy in four supergravity 
confining geometries; the Witten-QCD, the Maldacena-Nunez, the Klebanov-Tseytlin and the 
Klebanov-Strassler model. First we derived the PDE equations outlining the string prohle which is 
moving in the bulk of these models. Solving these PDE equations analytically, and then calculating 
the proper area between the event horizon on the prohle and the boundary can result in the free 
energy of accelerating particles. If one also finds the Unruh temperature of accelerating particles 
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in these confining supergravity backgrounds, then in the semi-classical limit one can find the 
entanglement entropy of pairs of particles moving on specific trajectories. 

We then studied the electric potentials in these models. By using the DB1 action we derived 
the critical electric field, and then by using a probe D3-brane, we derived the total potential and 
the distance between the quark and anti quark in each background. By solving the integrals 
numerically, we plotted the phase diagram of the total potential versus the distance between the 
quarks and anti quarks, x, in each geometry. As was expected we found three similar phases in 
each background. For the small electric fields or a no pair creation would be present. By increasing 
the electric field the Schwinger pair creation would happen by a tunneling process and with an 
exponential suppression. Finally by increasing the electric field the potential in all of these models 
becomes unstable catastrophically where the probability is no longer exponentially suppressed. 
The rate of this phase transition is faster for WQCD and KT, relative to KS and MN. We also 
presented the phase diagram of the Klebanov-Witten and AdS for comparing the results with the 
conformal cases. In the non-confining geometries the first phase does not exist. 

We then studied the phase diagram of the entanglement entropy of a strip in these confining 
geometries. We found the expected butterfly shape in the diagram of the entropy versus the 
length of the connected solution and a peak in the diagram of the length versus the minimal radial 
position in all of these backgrounds indicating the existence of a confinement phase. Again we 
found that the rate of changing phases is faster for WQCD and KT relative to KS and MN. It 
would be interesting to further study the possible relations between the Schwinger effect and the 
entanglement entropy. Finding such a relation could be very useful, as it can be used to measure 
the entanglement entropy in the lab with many applications in the condensed matter systems. 

In the final section we were interested to study the effects of an additional magnetic field on 
the rate of pair creation in each of these models. We found the critical electric field which makes 
the pair creation happening in the presence of the magnetic field and then by calculating the 
imaginary part of the Euler-Heisenberg effective Lagrangian, we studied the rate of pair creation 
for each model. The results was that similar to the Sakai-Sugimoto and deformed Sakai-Sugimoto 


models, 20 , the parallel magnetic field to the electric field increases the pair creation rate and the 
perpendicular magnetic field decreases the rate. This could be a universal feature of all confining 
geometries which can be demonstrated analytically in the future works. 
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